Exchange interaction effects in the thermodynamic properties of quantum dots 
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We study electron-electron interaction effects in the thermodynamic properties of quantum-dot 
systems. We obtain the direct and exchange contributions to the specific heat Cu in the self- 
consistent Hartree-Fock approximation at finite temperatures. An exchange-induced phase transi- 
tion is observed and the transition temperature is shown to be inversely proportional to the size of 
the system. The exchange contribution to Cv dominates over the direct and kinetic contributions 
in the intermediate regime of interaction strength (rs ~ 1). Furthermore, the electron-electron 
interaction modifies both the amplitude and the period of magnetic field induced oscillations in Cv 

PACS numbers: 73.21.La, 21.60.Jz 65.80.+n, 



I. INTRODUCTION 

Research on semiconductor quantum dots (QD)s and 
nanostructures have drawn considerable effort in recent 
yearsii In particular, the study of electron-electron in- 
teraction effects on the ground-state and excited-state 
of QDs has been a very active subject. A variety of 
methods have been used in such studies, ranging from 
the exact diagonalization of few electron systems^i^iii^iLi 
to sophisticated numerical schemes based on the den- 
sity functional theory, quantum Monte Carlo simula- 
tions and mean-field approximations!^ Among the last 
ones, the self-consistent Hartree-Fock (SCHF) approxi- 
mation has been successfully applied to QDs in a num- 
ber of ■ys^orks2ii24iiii24iiiiiiSii& which focused attention on 
calculations of the pair-correlation function^ and addi- 



tion spectr; 



.10.11 



and on configurations of the Wigner-like 



molecule in the strong interacting regimeiiSiii 

A less pursued track is the use of SCHF to study mag- 
netic and thermodynamic properties of semiconductor 
QPgj2^ji5ji6 Electron-electron interaction was shown to 
give an important contribution to thermodynamic prop- 
erties such as the magnetizationLi^ and the magnetic 
susceptibilityr-'-- Another quantity of experimental in- 
terest is the specific heat Cy, which has been studied 
in a number of works in both non-interactingiii^ and 
interactingS*15. QD systems. In Ref. 0, Dean and 
coworkers reported an interesting interaction-induced 
phase transition in parabolic QDs with iV ~ 6 electrons. 
This phase transition manifests itself as sharp drops in 
the specific heat as the temperature reaches a critical 
value. Nevertheless, a systematic study on how such a 
transition depends on the interaction coupling parame- 
ter of the dot, which measures the relative electron- 
electron interaction strength, remains to be performed. 

In this paper, we address the role of the exchange inter- 
action in the thermodynamic properties of non-parabolic 
QDs. Specifically, we study the kinetic, direct and 
exchange contribution to the specific heat in a finite- 
temperature Hartree-Fock approach^Siii In a previous 
work using this methodji^ we have shown that the ex- 



change interaction contribution is the dominant term in 
magnetic properties such as the zero-field susceptibility in 
the intermediate regime of interaction strength (r^ ^ 1). 
We find in this investigation that the exchange effects 
also play a dominant role on the specific heat properties. 
In particular, the exchange electron correlations domi- 
nate the finite-temperature phase transition and is the 
leading contribution to Cv for ~ 1. We also find that 
the transition temperature scales with the inverse of the 
dot size. As a consequence, this phase transition could 
in principle be experimentally observed for dots tens of 
nanometers across at an attainable temperature range. 

The paper is organized as follows. In Section ^1 we 
describe the system to be studied and its main features. 
The results for the specific heat and the discussion of 
main results are given in Section lllll as well as our closing 
remarks. 



II. THE MODEL 

We consider the problem of N interacting electrons 
confined in a 2D square quantum dot of size L and sub- 
jected to an external magnetic field B perpendicular to 
the electron system. To account for screening effects, the 
electron-electron interaction is modeled by an Yukawa- 
type potentials, and the model Hamiltonian reads as 



N 



N 



g2 ^-K\rn-r^^/ 



n<n' 



(1) 



where r„ indicates the position of the rtth electron. We 
consider low g-factor QDs, so that the Zeeman term can 
be safely disregarded. Above, k gives the effective inter- 
action range and is the background dielectric constant. 
For K = 0, there are no screening effects and the "bare" 
Coulomb interaction is recovered. 

The single-particle Hamiltonian h(r) is given by 



h{v) 



1 



2m* 



A(r) + u{v) , 



(2) 
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where m* is the electron effective mass and u(r) is 
the hard-wall confining potential. The vector poten- 
tial A is chosen in the symmetric gauge, namely, A = 
{—By/2, Bx/2,0). Hereafter, the magnetic field is ex- 
pressed in units of <i>/$Oi where $ — BA is the magnetic 
flux through the system area A and $o = he/ e is the 
quantum flux unit. 

A key parameter in our analysis is L/a*^, the QD 
length L in units of the effective Bohr radius — 
h'^er/rn*e^ which gives the relative strength of the e- 
e interaction as compared to the kinetic energy of the 
systemiiS For a square dot of side L, the potential en- 
ergy scales with while the kinetic energy scales with 
L~^. Therefore, as L is increased, the potential energy 
becomes increasingly more important. 

The standard dimensionless parameter that quantifies 
the ratio between the potential and kinetic energies of the 
system is the so-called Brueckner parameter , which in 
2D reads as = 
rs are related by r. 



A/iNT:[a 




Therefore, L/a*g and 
Furthermore, by 



choosing a square hard- wall confinement, one can easily 
tune Ts by changing the dot's lateral size L. 

The many-body ground-state energy is obtained in 
the finite-temperature SCHF approximation. The SCHF 
equations read asii*i^ 
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n, \ dr>*(r')«(r,r')</>,(r') 



E 



n, /drV*(r')«(r,r')0,(r)</),(r') 



-HF 



0.(r), (3) 



where the sums run over all HF orbitals. Here rii = 
{exp[(e^^ — /^)/fcBr] + l}~^ is the Fermi occupation num- 
ber of the ith HF orbital with corresponding wave func- 
tion </'i(r) and energy ef ^. As in the standard procedure, 
the chemical potential /i is determined by requiring that 
N ~ rii. We truncate the number of orbitals and take 
only the M > 2N lowest energy states into account. 
The SCHF ground-state energy is given by 



E 



HF 



E 



ni{(l)i\h\(j).. 



HF 



^^n^n■j(^{(|)^(f>j\v\(|),(|)j 



-{4'i4>j\v\4'j4>i 



(4) 



where the \(f>i) are the HF orbitals, self-consistent solu- 
tions of Eq. © and T^^ , Vj^^ and V^^ are the kinetic, 
direct and exchange contributions to the ground-state en- 
ergy respectively. We are interested in the intermediate 
interaction strength regime (r^ ^ 1). In this regime, the 
direct interaction term Vj^^ is the leading contribution 
to the Ef^, followed by T^'^ and V^^ respectively. 

The details on the matrix elements calculations and 
the Hartree-Fock method can be found on Refs. ^ and 

respectively. 

We calculate the specific heat C„ of the system in this 
SCHF approximation, namely: 



df)y 



(5) 



where U is the internal energy and T is the system tem- 
perature. The first-order exchange and Hartree contri- 
butions to Ct, can be accounted for by approximating 
U ~ E^^ so that for E^^ , given by Eq. 0] there are 
kinetic (C|^'"), direct (C^) and exchange (C^) contribu- 
tions to Cy. 



III. RESULTS AND CONCLUSIONS 

We analyze the behavior of the specific heat Cy as a 
function of relevant parameters of the system, i.e. the 
temperature T, the interaction strength L/a*^, the mag- 
netic field <&/<&0j the potential range and the num- 
ber of electrons N . The results given in this section 
are for the Coulomb case (k = 0) with iV = 10 elec- 
trons in the dot. We should mention that QDs with up 
to iV = 20 electrons and with the screened interaction 
(k 7^ 0) were also considered and the same overall qual- 
itative features were observed. In the remaining of this 
section, energy and temperature are given in units of the 
typical scales for the system, namely E^ ~ / {m* L'^) 
and Tl = E^/kB respectively. 




FIG. 1: (color online) (a): Ground-state energy -Ehf as func- 
tion of temperature for A'' = 10 electrons and rs = 0.89. The 
kinetic (b), exchange (c) and direct (d) contributions to iSjiF 
are also plotted. 



The ground-state energy E^^ increases with temper- 
ature as shown in Fig. ^ for A = 10 and rg — 0.89. 
Nevertheless, this increase is not smooth and a sudden 
change in slope is observed at a certain temperature T* . 
By analyzing the kinetic energy T^^, the direct term 
Vj^^ and the exchange contribution given by Eq. 
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FIG. 2: (color online) Specific heat as function of temperature 
for the noninteracting (solid line) and interacting cases. In 
the latter, the interaction strength values are Va = 1.07 (up 
triangles), 2.14 (circles) and 3.30 (diamonds). The insets show 
typical charge distributions in the dot before and after the 
transition. 

we see that both Vj^^ and V;^^^ decrease with T, 
as expected. Furthermore, we observe clearly that the 
change in slope is a feature due to the exchange interac- 
tion, since neither T^^ or Vj^^ display cusps at T = T* 
(Figs. IBj-d). 

The abrupt change in slope in the energy causes a dis- 
continuity in the specific heat CyiT) at T = T*, as seen 
in Fig. 121 A sharp drop develops for a wide range of val- 
ues of the interaction strength parameter (for simplicity, 
the noninteracting case is referred to as "L/a*g = 0" or 
"r^ = 0" . The noninteracting curve does not display any 
sharp drops). 

Discontinuities in the specific heat are usually regarded 
as signatures of phase transitions.^^ In fact, such tran- 
sitions are accompanied by a charge reordering in the 
ground state distribution, as shown in the insets of Fig. 
121 The lack of the C4 rotational symmetry in the charge 
distributions is a consequence of the nonlinear coupling 
of the original orbitals in the Eq. (PJ and is a peculiar- 
ity of the Hartree-Fock approximation. Nonetheless, a 
clear charge rearrangement is verified as the system un- 
dergoes the phase transition. Most strikingly is the fact 
that this is an exchange-induced phase transition and it 
is a direct consequence of exchange effects between the 
electrons in the dot. Such transitions were reported in 
previous studiesi^ as related to phase transitions in the 
ground-state charge distribution inside the dot. 

The reduced transition temperature varies both as a 
function of the number of electrons and the relative inter- 
action strength. For N = 20, up to three transitions are 
observed in the temperature range 2 < T/Tl < 20 (not 
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FIG. 3: Reduced transition temperature T* /Tl as function 
of (see text for details). Inset: T* /Tl increases linearly 
with the interaction strength L/ag. 



shown). In Fig. O the dependence of the critical temper- 
ature T* /T]^ with (which is proportional to the den- 
sity N/L^) is shown. We also depict T* /T^ as a function 
of L/a*g in the inset of Fig. Oand a roughly linear depen- 
dence is observed. These results imply that the transition 
temperature T* scales with since Tl oc L~^. The 
data are very well described by a linear fit (T* cx L^^) 
as shown in the figure. 

Such a scaling behavior allows us to estimate the values 
of the critical temperature T* for typical dot sizes. For 
GaAs QDs with L ^ 50 nm one obtains T* ~ 11 K 
and it decreases with 1/L for larger dots. The transition 
temperature further decreases when screening effects are 
taken into account. For the potential range of = 
L/10 and L - 50 nm, we obtain T* ~ 2 K. 

We have also investigated the specific heat dependence 
with both the interaction strength parameter and the 
magnetic field for a fixed temperature. The kinetic, di- 
rect, and exchange contributions to Cy a,t T = 3Tl are 
shown on Fig. 01 as a function of . As the relative 
interaction strength parameter increases, the exchange 
contribution rises fast and becomes the leading con- 
tribution to Ci, for Vs 1. This is a surprising result 
since V^f^^ is smaller than Vj^^ and T^^ by a factor of 
2 — 3 (see Fig. However, the effect of temperature 
in the variation of V^^ is stronger and > C^™, 
for Ts ^ 1. The direct term, on the other hand, gives a 
smaller negative contribution which cancels out the pos- 
itive C^'" (which is dominant for ^ 1). A peak in C„ 
appears at the value of Ts for which T* = 3Tl and ap- 
proaches zero for higher values of r^, since temperatures 
changes do not sensibly affect the ground-state energy in 
the strongly interacting regime. 

The specific heat oscillates as a function of the mag- 



4 




12 3 

r 

s 



FIG. 4: (color online) Specific heat as function of (circles). 
Also plotted are the kinetic (squares), direct (up triangles) 
and exchange (down triangles) contributions. Inset: C„ as 
function of magnetic flux ^/^o for the noninteracting (solid 
line) and interacting cases with — 0.53 (triangles), 0.89 
(half- filled circles) and 1.78 (filled circles). Curves are off-set 
for clarity. 

netic field with increasing amplitude, as seen in the inset 
of Fig. 01 The interaction influences both the period 
and the amplitude of C^iB). The noninteracting curve 
displays an oscillatory pattern with both high and low 
harmonics. For higher values of the interaction strength 
parameter, the higher harmonics are suppressed and an 
oscillation period is defined more clearly . For even higher 
values of rs, the oscillation period decreases. A similar 
behavior was seen in the magnetization and magnetic 
susceptibility of QD systcms^^ and it is related to an ef- 



fective increase in the chemical potential as increases. 

In summary, we have investigated interaction effects 
in the thermodynamic properties of QDs. The exchange 
interaction plays a relevant role on the specific heat fea- 
tures and is the leading contribution for dots in the 
~ 1 range. The exchange induced finite-temperature 
phase transition, studied in previous works, ^- has been 
shown to depend on the interaction strength parameter, 
or equivalently, on the system size L. The transition tem- 
perature T* decreases as L increases and we estimate that 
the transition regime could be experimentally accessible 
for sufficiently small dots. Furthermore, the specific heat 
oscillates with the magnetic field and both period and 
amplitude of such oscillations strongly depend on inter- 
action effects. 

Experiments to verify our findings using single 
quantum-dots are likely too demanding. The specific 
heat has been measured in multi-layer 2D electron gases 
in the Landau regime with heat-pulse^ and steady-state 
ac-temperature calorimetrySl techniques, with resolu- 
tions in Cy still much lower than the required to test 
our results. One possible way to overcome such difficulty 
is to perform experiments in ensembles of nearly identical 
dots in a multi-layer configuration so that the contribu- 
tion from single dots is amplified. This is, nonetheless, 
an experimentally challenging task which would bring a 
new understanding to the many-body effects in the ther- 
modynamics of such small-scale devices. 
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